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Abstract 

Light-front QCD is studied by the method of flow equations. Dynamical gluon 
mass is generated, which evolves with the cut-off according to renormalization 
group equation. Eliminating by flow equations the quark gluon coupling with the 
dynamical gluon mode, one obtains an effective interaction between quark and 
antiquark which exibits the Coulomb and confining singularities. The scale, which 
regulates the light-front IR singularities in the gluon sector, defines the string 
tension of confining interaction. The mechanism of confinement in the light-front 
formalism is suggested, based on the singular nature of the light-front gauge. 
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1 Introduction 



Quantum Cromodynamics (QCD) is a widely accepted theory of strong interactions. 
This wide-spread acceptance is based on the success of Feynman rules of perturbative 
covariant calculations, which provided convincing agreement between perturbative QCD 
and experiment. However there is a gap between the perturbative behavior of QCD and 
its low-energy limit, where perturbation theory breaks down and the physical observables 
such as mass spectrum and decay width are predicted based on phenomenology. In the 
strong coupling regime non-perturbative methods are required. The non-perturbative 
solution of a bound state problem can be obtained directly only by using the Hamiltonian 
formalism. 

It is crucial to a successful non-perturbative solution that it exposes the three impor- 
tant long range properties of QCD: confinement, spontaneous chiral symmetry breaking, 
and the topological structure. We show in this work, that the non-perturbative method 
of flow equations, when applied to the QCD Hamiltonian in the light-front quantization, 
provides some understanding of the physics of confinement in the Scrodinger picture. 

There were several successful attempts to reveal the mechanism of confinement in the 
Schrodinger picture, using a special ansatz for a vacuum wave functional and integrating 
over all possible gauge configurations [|IJ. We believe that the light-front quantization 
provides an alternative formalism, where it is possible to isolate the degrees of freedom 
that are responsible for the long-range properties of QCD. In the light-front quantized 
QCD the topological structure is carried by the zero mode of A + instead of (nontriv- 
ial) gluon vacuum configurations as in other gauges. 

The confinement mechanism in the light-front QCD was suggested several years ago 
PI, based on the fact that QCD 3+ i already has a confining interaction term in the 
light-front Hamiltonian, the instantaneous four Fermion interaction, which is the con- 
fining interaction in QCD 1+ i. The authors argue, that in QCD the second order quark 
glue interaction, which appears through similarity renormalization, does not cancel the 
instantaneous interaction as it does in perturbation theory. The singular part of the un- 
canceled instantaneous interaction (~ l/q +2 ), produces a logarithmic potential of the 
form 



where a is equal 1 for the radial tensor e r along the z-axis and it equals 2 when f*is purely 
transverse; a> has dimension of energy and sets up a scale for the potential. This potential 
is confining. It is boost invariant, but it is not rotational symmetric that is confusing. 
It was also pointed in that the instantaneous interaction of the off diagonal currents 
is modified by the topological properties of the theory, and confinement is destroyed for 
these currents in QCDi + \ and also in QCD 3+ i. 

The basic idea of our calculations is to use the 'singular' nature of the light-front 
gauge in conjunction with non-perturbative renormalization provided by flow equations. 
The main conceptual complication to study renormalization group (RG) flow of the ef- 
fective QCD Hamiltonian is the lack of a well defined initial condition. We study the 
RG equation for an effective gluon mass, given some parameter mass fi for the gluon 
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mass in renormalization point. This mass p is used as a parameter and is taken to zero 
at the end of calculations. Apart from the known perturbative gluon correction in the 
second order, we obtain also the term which is a function of the ultraviolet cut-off A 
and the mass parameter p. We associate the second term with 'non-perturbative' gluon 
mass correction. In the light-front quantization the gluon non-abelian gauge interactions 
produce severe infrared divergences in the effective gluon mass. The subtle point is, that 
we use an additional scale m, m C A, as suggested by Zhang and Harindranath in |J, 
to regulate these divergences. One ends up with an effective 'non-perturbative' gluon 
mass, which is the function of the parameters u and p and the cut-off A, /i|r P (A; u, p). 
In the limit p — > 0, only the non-abelian part of the effective gluon mass, regulated by 
the scale u, survives. This is the crucial difference between QED and QCD. In QED the 
perturbative gluon mass correction can be removed by perturbative renormalization, i.e. 
renormalized photon stayes massless. In QCD by absorbing the leading cut-off depen- 
dence in the second order mass counterterm, we are still left with non-perturbative mass 
correction. It turns out that this 'dressing' of an effective gluon playes an important role 
in the effective interaction between quarks. 

We simulate an effective interaction in QCD between probe static quark and antiquark 
as an exchange of 'non-perturbative' gluon (gluon flux) with a nonzero effective mass 
//jvp(A; u, p) (effective energy), which evolves with the cut-off A according to RG equation. 
Eliminating by flow equations the quark gluon coupling with an effective gluon mass, 
one obtains the quark-antiquark potential which includes two pieces: the short-range 
part describes the perturbative one-gluon exchange and is analog of the perturbative 
interaction in QED |5|]; the long-range part arises due to the non-perturbative gluon 
'dressing', i.e. due to the dependence of gluon effective mass on the cut-off. For the 
vanishing mass parameter, p — > 0, an effective gg-potential is given by a sum of Coulomb 
and linear rising confining interactions 

rv 

V(f) = -C f — + a-r (2) 

r 

The parameter u sets up a scale for the long-range part of interaction: it defines the 
string tension of the confining term, a ~ u 2 . Though the calculations are performed in 
the light-front frame, the resulting effective interaction manifests rotational symmetry. 

The article is organized in the following way: the first section introduces the problem 
and sets up a scheme; in the second section a gluon gap equation is obtained and solved 
for an effective gluon mass; we obtain an effective potential between quark and antiquark 
in the third section. 

2 Flow equations in QCD 

Flow equations were discussed in great detail in application to QED in the previous 
work ||. Here we point out the differnce between QED and QCD. For simplicity we 
consider only the abelian part of QCD Hamiltonian. Flow equations for the Hamiltonian 
H = H d + (H- H d ) read i 

^ = fo(0.*(01 
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V (l) = [H d (l),H(l)}, 



(3) 



where eta is the generator of unitary transformation, which eliminates the particle num- 
ber changing part of the Hamiltonian, (H — H^); the H d includes all particle num- 
ber conserving terms; I is the flow parameter, which changes from I = correspond- 
ing to the initial canonical Hamiltonian to / — > oo with block-diagonal Hamiltonian 
H(l -f oo) = H d {l -> oo). 

It is always possible to divide the complete Fock space (particle number space) into 
two arbitrary subspaces, P and Q space. The Hamiltonian matrix reads 

/ PHP PHQ \ 

\ QHP QHQ J ' w 

where P and Q = 1 — P are projection operators. For the (abelian) QCD the content of 
sectors is given 

PWi = \g) 

W> = \qq), (5) 

with simbols g and q standing for gluon and quark, respectively. Therefore the matrix 
elements of PHQ descibe quark gluon coupling, PHP stands for gluon effective energy, 
and QHQ descibes qq effective interaction. When the Hamiltonian matrix is subject to 
the unitary transformation Eq. (^), the secor Hamiltonians become the functions of the 
flow parameter /. 

Suppose we know approximately the eigenstates of the sector Hamiltonians PH(l)P 
and QH(l)Q and their eigenvalues E p (l) and E q {f). The indeces p and q run over all states 
in the P and Q space, respectively. Suppose further, that this basis is /-independent, i.e. 
we assume, that the off-diagonal matrix elements h pp i and h qq > of PHP and QHQ are 
small. For the particle number conserving sector we keep all the terms in flow equation, 
while for the particle number changing sector we neglect the small off-diagonal matrix 
elements h pp i and h qq ' and take into account only the diagonal matrix elements E p and 
Eg on the right-hand side of flow equation. Flow equations for the matrix elements of 
the particle number conserving and particle number changing sectors read, respectively, 



dh ppl (l) 

dl 
dh pq (l) 

dl 



-(E p (l)-E g (l)) Vpq (l), (6) 



and the analagous equation for h qq i. Here the generator is chosen in a more general, than 
Eq. (0), form 

z pq (l) = l{E p {l)-E q {l)f , (7) 
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where f(z) is the similarity function with the properties 

/(0) = 1 

/(z-oo) = 0. (8) 

We take into account in the similarity factor the dependence of the energies on the 
flow parameter. This is the crucial difference between QED and QCD. We show in the 
next section that E p (l) playes the role of the effective energy (effective mass) of gluon. 
In QED this dependence can be removed by perturbative renormalization, so that one 
works in terms of renormalized energy (mass) operators which are fitted to the physical 
values. In QCD only the perturbative energy (mass) correction can be absorbed by the 
counterterm. The non-perturbative energy correction, which is left, shows how gluons 
(and quarks) are getting 'dressed' from bare to constituent degrees of freedom. 
The solution for the particle number changing part reads 

WO = Wo)/(W0) > ( 9 ) 

which shows, that only matrix elements in the band \E p (l) — E q (l)\ < 1/y/l = A survive. 
In the third section we consider possible choices for the similarity function. For example, 
for the similarity function 

f(z) = exp(-v^) , (10) 

the particle number changing part decays exponentially as I — > oo (or A — ► 0). In the 
case of degenerate eigenvalues of initial Hamiltonian, E p (0) = E g (0), the paritcle number 
changing part still decays, but algebraically 

exp(-v^) ~ (A//2)" 2A/A , (11) 

due to non-perturbative dependence of energy on flow parameter 5E NP (X) ~ /iln(A 2 //i 2 ), 
where jl is some energy scale (see the second section). For the particle number conserving 
sector one has 

dh pp ,(l) _ ^ / dh pq (l) 1 1 dh qp ,(l) \ 

dl ~ dl E p (l)-E q (l) qp ' [)+pq[) E p ,(l)-E q (l) dl J' [ } 

and analagously for the Q space. Here h pq (l) is given by Eq. (|S|). When the sectors are 
assigned as in Eq. @, the equation for the diagonal matrix elements in P space, p = p', 

^--E^y^l(MOMO), d3) 

provides (after integrating over the flow parameter) the gap equation for an effective 
gluon mass. The equation in Q-space 

dhggjl) \ ^ / dhgpjl) 1 1 dh pq ,(l) \ 

dl ' " dl E q (l)-E p (l) pq ' [)+qp[) E q/ (l)-E p (l) dl )' K ) 

defines an effective qq interaction. The ultimate aim is to solve these equations selfcon- 
sistently. In the next two sections these equations are solved analytically doing some 
approximations. 
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3 Gluon gap equation 

Integrating flow equations over the flow parameter in one-body sector gives gap equations 
for the effective energies of quark and gluon, Eq. Provided the connection between 
light-front energy and mass is given p~ = CO f or quark and q~ = 9±+ q + f° r gluon, 

flow equations for quark and gluon effective masses are 

dm 2 (i) rdkid 2 k 1± 2m i df'(D 3 (iy,i)o(ki)e(kt) 



dl ^ ^ 'J 16tt 3 * qK * J D 3 (l) dl kf kt 

x u^D^hWifa + m(l))Yu(p)5^(p -fa- k 2 ) , (15) 

and 

V(0 x T (TaTb , f dkfd 2 fa ± 1 dp(D 2 (iyj) Q(kf)Q(kt) 
-dT 9 ^ ~ " Tr(TT) i 16vr3 9 ^ l) D^) dl k? W~ 

x Tr ( 7 "(ft + m(l))Y(- fa + m{l))) 5 {3 \q -fa- fa) (16) 
rdkt<Pk 1± , n , 1 dp(Di(l);l)&(kf)&(kt) 



x 



2 (f aC f I] ^3 9g(l)- Di ^ ^ k + k + 

TW(q, -fa, -fa)D acT ,{fa)Y v ^' {-q, fa, fa)D p r p (-fa)6®(q -fa- fa) 



where in the last term the factor | is a symmetry factor for two-boson states. Note, 

that no correction arises to the term ^ (and which is protected by the kinemat- 
ical symmetries; the total transverse momentum does not appear in a boost invariant 
expression. 

The gluon couples to the quark- ant i- quark pairs and pairs of gluons while the quark 
couples only to the quark-gluon pairs. Here the energy differences are 

Di = q~ — ki — k^ 
D 2 = q~ — k^ — k% 

D 3 = p~ - k^ - k 2 . (17) 

One should not be confused, that the momenta in different loops are denoted by the 
same letters, fa and fa. The trigluon vertex is (— g g )T pup , with || 

T^(p, q, k) = (p- qfgr + (q- k) p g pv + (k - pfg^ . (18) 

Generally, the quark-gluon coupling, g q {l), and the trigluon coupling, g g (l), are different 
from each other (functions of three vertex momenta) for a nonzero flow parameter. The 
energy differences depend on the flow parameter I through the masses, i.e. D 3 (l) = 
[p]_ + m 2 (l))/p + — [k\ ± + fj, 2 (l))/ki — (k% ± + m 2 (l))/k2. The polarization sum (in the 
light-front gauge) is given || 

DAk) = E e^KW = ~9r + ^ K l^ K , (19) 

A=l,2 K 
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where k ■ e = 77 • e = 0, and the light-front vector 77 is defined as rj ■ k = k + . The sum over 
helicities for the Dirac spinors is 

u(p, a)u(p, a) = $ + m 

<T=±l/2 

v (p, a)v(p, a) = p 1 — m. (20) 

(T=±l/2 

We introduce 

Qf(X) = —q + D 1 (X) 
Ql(X) = -q + D 2 (X) 

Ql(X) = -p + D 3 (X). (21) 
In the light-front frame Eq. (pl|) and Eq. fllTD read 

dm 2 (X) fi dx r d 2 k ± 1 df 2 (Q 2 (X);X) 



and 



dX 1 Jo x(l-x)Jo lQTt^ qy J Ql(X) dX 

x fe(-^- + 4) + 2m 2 (A)-^) , (22) 
V 1 — x x z 1 — x 



dfi 2 (X) 0T fdx r d 2 k ± 2 1 d/ 2 (Q|(A);A) 



dA 1 J Jo x(l-x)Jo i6tt w Q\{X) dX 

( k 2 ± + m 2 (X) \ 
\ x(l-x) L J 

^ f dx r d 2 k ± 2 1 # 2 (Q 2 (A);A) 



x(l-x)7o 16vr 3ysV ; Q?(A) dX 
x (1 + 1 + 7-^-^)1 , (23) 



x 2 (1 — 



where 



2(1 — xj 

Q 2 (A) = fci + m2(A) + fci + /i2(A) -m 2 (A), (24) 
a; 1 — x 

and we used the connection between the flow parameter, /, and the ultraviolet cut-off, 
A, as I — 1/A 2 . Here Casimir operators in fundamental and adjoint representations are, 
respectively, C f = T a T a = (N 2 - l)/2N e and C a 5 ab = f acd f bcd = N c 6 abl N c is the number 
of colors (i.e., N c = 3); and Tf5 a b = Ti(T a T b ) = \5 ab . This system of equations, Eq. (f22j) 
and Eq. (|23|) , was considered also in [0. 
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Generally, it is very difficult to solve both equations self-consistently. One of the 
reasons is that the equations involve running coupling constants, g q (X) and g g (\), which 
depend on A in accordance with renormalization group equations. (RG equation for the 
couplings can be obtained from the flow equations for the quark-gluon vertex in the 
third order and for the trigluon vertex in the forth order, respectively). Also initial 
conditions for these equations are not known. In the leading order one can decouple 
the gap equations for quark and gluon effective masses and the RG equation for the 
coupling constants if we neglect all dependences on the cut-off in the right-hand side of 
the corresponding flow equations. In order to take into account these dependences we 
need to consider higher orders. 

Below we study the gluon gap equation, that generates an effective gluon mass which 
depends on the cut-off A. In the next section we show that the exchange by a gluon 
with the cut-off dependent mass leads to a confining potential between static quark and 
antiquark at large distances. 

We integrate the flow equation for the gluon energy, Eq. ( ^3|) , neglecting cut-off 
dependence of masses and coupling on the r.h.s. The constant of integration is assumed 
to be a 'physical' mass (for derivation see JF|). We take for an effective gluon mass some 
value ft. Gluon gap equation for the abelian part (quark loop) reads 

//(A) = ^ + 2g%N f J o 1 dxJ o C °^f(QlW/^ 



/i 2 (A)(2a; 2 -2x + l) 2m 2 



\k]_ + m 2 — x(l — x)n 2 (\) k]_ + m 2 — x(l — x)/x 2 (A) 

+ <- 2 + ^)))- (25 > 

where we rescaled the cut-off A — > X 2 /q + . Similarity function plays the role of UV cut-off 
in the loop integral. It also regulates the light-front IR divergences (for m^O) 

k ± ma x = x{l - x){\ 2 + ^ 2 (A)) - m 2 



2 2 

m m 



V + ^(A) 2l21 -]^(A)' (26) 



One has 

,2 , 9 2 T f N f 



^(A) = jl 2 + 9 \ f dx U 2 (A)(2x 2 - 2x + 1) In 

07T J%min \ 



x(l — x)X" 



m 2 — x(l — x)jji 2 {\) 



+ 2m 2 In 



x(l — x)A 



2 



m? — x(l — x)n 2 (X) 
A 2 + /x 2 (A) 



+ jj(A 2 + /AA)) 



- 2m 2 (ln ■ (27) 

m z J 

When the renormalization point is taken at q 2 = and jl 2 = the gap equation Eq. 
(|25D (and Eq. (|^)) is reduced to the perturbative case. The perturbation correction is 
given 

^„(*) = «£. (28) 
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Non-perturbative solution of the integral equation Eq. fl2"7|) can be obtained numeri- 
cally. Instead, we solve Eq. (|27f) itteratively. In the leading order /U 2 (A) = ft 2 . The next 
order reads 



f 2 m 



■Sit' 



dx [ft 2 (2x 2 -2x + l)ln 



x(l — x)X 2 



m 2 — x{l — x)/i 2 



+ 2m 2 In 
- 2m 2 (In 



x(l — x)\ 2 



m 2 — x(l — x)jji 2 

nV 



+ ^(A 2 +/i 2 ) 



X 2 + ft 2 



(29) 



This is equivalent to take the renormalization point q 2 = fi 2 . Provided that m<C/i<C A, 
Eq. (|29| ) is reduced 



~, 2 + <gW f A 2 



4vr 2 



~.2 



/i' + S/ip T (X) + 5/j? NP (X;fl) . 



(30) 



We add the finite part (independent on the cut-off A) to ft, the result is denoted as ft'. 
We associate all the terms which depend on A but do not depend on ft with perturbative 
correction, denoted as <5/Xp T (A). The rest, except constant ft,', gives non-perturbative mass 
correction Sft 2 NP (X, ft). The perturbative term is given in Eq. (p8|) . 

We consider the non-abelian part (gluon loop), Eq. (p3[). For simplicity we consider 
exchange with massless gluon; the external gluon is put on mass-shell q 2 = /i 2 (A), i.e. 
the energy denominator is 



or (A) 



x(l — x) 

Gap equation (non-abelian part) reads 

r °° d 2 k 



aAa) • 



M 2 (A) 



2g 2 C a f 1 
Jo 



dx 



16tt 3 



x 



1 + 



/j, 2 (X)x(l — x) 
k 2 ± -x(l-x)fi 2 (X) 



f(Q\ 2 w/y 

1 x 2 ~^ 



l-x] 



The range of integration is defined by the similarity function 



k 2 



(A 2 + fi 2 (X))x(l-x) 
< x < 1. 



(31) 



(32) 



(33) 



The IR singularities, from x — > and not regulated by the UV cut-off for kj 

since the massless gluon is exchanged. We regulate it by principle value prescription | 



1 

x 



1 1 

+ 



2 \x + ie x x — ie 2 



(34) 
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where e x is dimensionless and is boost invariant (i.e. e x = e/q + ). The result of integration 
of Eq. (^) reads 

/^(A)=A' + ^((A' + AA))(h.i- jij+AAJb^Ch,!- jl)) . (35) 
The perturbative correction, when q 2 = and fi = 0, is given 

The itterative solution reads 

= /? /2 + ^ T (A) + ^ P (A;/i), (37) 

the finite constant is denoted as //", and /i" ~ /L 

We combine quark Eq. (|30D and gluon Eq. ( p7|) loops to give the leading order result 

= p 2 + 5[f PT (\) + 5fi 2 NP (X;fi) , (38) 
where the perturbative term 

^pt(A) = £~ 2 X 2 (C a (\n I - H) + r /JV/ i) , (39) 

reproduces the known result for the perturbative mass correction (without instantaneos 
graphs) 0. Here we assumed for the constant term p! + \i" ~ /i 2 . 

The (non-perturbative) gluon mass correction S/j,% P (\,p,) contains logarithmic IR 
divergence. Even by adding instantaneous graphs one can not eliminate severe IR diver- 
gences that appear in the gluon sector []||. One may introduce a non-zero mass jx for 
gluon in the intermediate state to regulate these divergences (Appendix A). But this will 
cause a mass singularity from massless gluon parameter \x — > for the on-shell gluon 
q 2 = f|. If gluon mass parameter in intermediate state, introduced as IR regulator, 
is the same as the mass of in- and out-going gluon (i.e. q 2 = ft 2 ), then an effective 

1 There are also instantaneous diagrams, which arise from the normal-ordering Hamiltonian, and 
in principle, they should accompany the generated (dynamical) terms in flow equations. Following the 
rules of light-cone perturbation theory S , one can take into account instantaneous graphs by replacing 
intermediate momenta in the nominator of dynamical diagrams as k = (k + , J2 in k~ — ^2' interm k~, k±). 
For example, for the quark effective mass, Eq. (|l^), in order to add instantaneous gluon exchange one 
should make the change k\^ — > k\^ = — fc^; to add instantaneous quark exchange k2 fl — > k^ = 
Pp. — fci/j- Analagously in gluon sector, Eq. ([L7l) . 
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(36) 



(non-perturbative) gluon mass vanishes as fx — > (Appendix A). These problems can be 
avoided if we introduce a mass scale u, as suggested by Zhang and Harindranath ||, for 
the minimum cut-off for transverse momentum, k±. This is equivalent to the integration 
of flow equitions for an effective gluon mass over the flow parameter in finite limits {u; A} 
[0. The integration of Eq. (the non-abelian part) gives 

/i 2 (A) = Lc 2 (u) + 2g 2 C a dx ?^(/ a (0?;A)-AG?;«)) 



x 



JO 16/T 

fi 2 x{l — x) fi 2 



kj_ — x(l — x)fi 2 + fi 2 k\ — x(l — x)fi 2 + jl 2 / 
1 1 
x 2 (1 — x) 



x ( 1 + -F + 7T— • ( 4 °) 



where the renormalization point is q — fx , and 



P 4- n 2 

Ql = ^ ±J ^-f- (41) 

x(l — x) 



The similarity function restricts the transverse momentum to 

k±max = (A 2 + fl 2 )x(l - X) - fi 2 



k±min = (U 2 + fL 2 )x(l - x) - fi 2 , (42) 



and the x range of integration, provided u -C A, is 



<x<l--^-. (43) 



■u 2 + fl 2 U 2 + fi 2 



Integration of Eq. (f40|) gives 
„»(A) - „»(„) = (A' h, J + h, J - A) + (A 2 - f 2 - ] , (44) 



where /i 2 <C m 2 . The value /i 2 (w) can be found from a renormalization condition for 
the 'physical' gluon mass 0. The following condition for fitting fi 2 (u) is assumed: the 
effective Hamiltonian at the scale u has bosonic eigenstates with eigenvalues of the form 



q =^^, i.e. 



where /i denotes the 'physical' gluon mass; 77 ^ is the first order generator which elim- 
inates the quark gluon (trigluon) coupling constant H^; \q) denotes a single effective 
gluon state with momentum q + and q±, (q'\q) = 167r 3 q + 5^(q' — q). In fact, the initial 
gap equations, Eq. (|25| ) and Eq. (0), were obtained using this renormalization condition, 
Eq. (|45|). In the order 0(g 2 ) one obtains 

~2 2/x\ /_ 2 , A 2 , U 2 U 2 5 2 U 2 11 A 

/i =MA)- 1 ^U in^(-^ + in^--) + A(in^--)l • ( 6) 
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So, 

Au) = ^ + 9 4^(^^t 2 (-^ + \nt 2 -h+u\\n%~ 1 ^)) . (47) 
47r 2 \ /i 2 n 2 n 2 12 /i 2 12 J 

Therefore, the effective gluon mass in the interacting Hamiltonian at the scale A is given 

" 2 w^ 2 +#(^'V ( -^ +1 V-^ )+A2(1 V-^) • (48) 

where parameter u 2 plays the role of IR regulator. This equation can be compared to 
Eq. 

To be consistent, we integrate also the quark loop over the flow parameter in finite 
limits, though in this case the IR singularities are regulated by the nonzero quark mass, 
m ^ 0. In full analogy with nonabelian case, one has from Eq. ( p3[ ) (the abelian part) 

^(A)-^) = ^(^ 2 1.4 + m2| V + ^ (A2 - t ' 2) ) ' (49) 

Therefore 

2/^ ~2 9 2T f N ff 1 ~2, A 2 2l A 2 1 2 \ 

/i 2 (A) = f + V-J-L (^ 2 In - + m 2 In - + -A 2 J , (50) 

that one can compare to Eq. (^). We combine the abelian, Eq. (j50D, and nonabelian, 
Eq. (p|), terms 

= ^ + ^ 2 fa(ln^-li )+r/iV l X 



4vr 2 V P? 12 3 / 
= /i 2 + <5/ip T (A) + 5fj 2 NP (X,fi, u) . (51) 



The sum of constant terms in Eq. (|50|) and Eq. (|48|) is of order of the mass parameter 
jl 2 . Here the perturbative term is 



Vpt(A) = ^A 2 (C^ln - - -) + TfNf- j . (52) 

Note, that for «> jl the nonperturbative mass correction in Eq. ( |5"I| ) changes the sign as 
compared with nonregulated mass correction Eq. (0), since the intermediate gluon with 
nonzero mass contributes negative term (see Appendix A or Eq. (flUD ). In fact, QCD with 
nonzero gluon mass resembles Yukawa theory. It was shown for Yukawa theory, that the 
leading correction proportional to A 2 and the logarithmic correction ~ lnA 2 //i 2 appear 
with opposite signs 

We introduce the mass counterterm to renormalize the Hamiltonian perturbatively 
in the second order 

m CT = ~ <^A*pt(A = A — > oo) , (53) 
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that removes the leading cut-off dependence. The rest is the (nonperturbative) effective 
gluon mass 

A 2 

VnpW = jl 2 - &(jl,u) In — , (54) 

where we introduced 

^)--^i^4 +, 4-^ +T ' N 4 + P) ■ < 55 * 

The limit of the <r(/i, u) as jl — > is finite, and is equal 

a = lim a(jl, u) = -^u 2 , (56) 

that plays the role of the string tension between quark and antiquark (see next section). 
One should not be confused that the string tension is proportional to the coupling con- 
stant. In Appendix A we show, that by the proper regularization of IR divergences the 
string tension has a pure non-perturbative form. In Eq. fl55|) we can remove the coupling 
constant by rescaling u 2 — > u 2 /g 2 , where g is the renormalized coupling constant. 

In Eq. fl54|) for A <C jl the effective (nonperturbative) gluon mass equals the 'physical' 
gluon mass, mass parameter jl. For A>/i there is the nonperturbative correction to the 
mass jl, given by the second term, which describes 'dressing' of gluon. 



4 Confinement 

Eliminating the quark gluon coupling one obtains an effective interaction between quark 
and antiquark, Eq. (]H|). The effective interaction between electron and positron in the 
light-front gauge was obtained in the previous workQ || 

V eS = -4n 2 a em (rr)B^, (57) 
where the current-current term in exchange channel is given 

(yv) = W± A'lhMpi, AQ) (vfa, A 2 ) 7 M?4 ^ 2 )) ^ (5g) 

^xx'{l — x){\ — x') 

where x = pf / P + is the longitudinal momentum fraction. The energy transfers along 
electron and positron lines are given, respectively 

Di = PT - p'i - q~ 

D 2 = p' 2 ~ P2 - <T , ( 59 ) 



2 The difference in the prefactor between ||] and Eq. ( |57| ) is x (j^j) , which comes from the light-front 
normalization in the bound state integral equation ~ J j^t- 
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where q = p\ — p[ = (q + , q±) is the photon momentum. The energy differences Eq. (^9|) 
are related to the four-momentum transfers 

Ql = -(Pi-Pi? = -q + D l 

Ql = -(P2-P2) 2 = -q + D 2 . (60) 
The following combinations are useful 

Q 2 = (Ql + QD/2 

SQ 2 = {Ql-QD/2. (61) 
We introduce also the energy differences in t-channel 

D[ = p^ +p 2 - (pi +p 2 y 

D' 2 = pT+p' 2 --(p[+p' 2 r, (62) 

and will use them later. They are related to the invariant mass-squares of the initial and 
final states 

Ml = (pi+P2) 2 =p + I>i 

M 2 2 = (p[+P 2 ) 2 =P + D' 2 , (63) 

also 

M 2 = (Ml + Ml) i '2 
5M 2 = (Ml -Ml) j2. (64) 

The tensor part B^ u of the effective interaction Eq. (B7I) includes two terms 

B, v = B*? + B; n u st . (65) 

The first one is generated by flow equations in the second order of perturbation theory 



*r = («) (§§+|§ ) . ( 66 ) 



and describes the dynamical photon exchange between electron and positron in s-channel. 
Here the polarization sum D^ u is given in Eq. (|T9"D; and the energy denominators are 
given in Eq. (0). The G-factor is defined 

Q^fdX ^^ fiQim, (67) 

where in order to preserve the boost invariance the cut-off is given in units of longitudinal 
exchange momentum, i.e. X 2 /q + . The meaning of the B-factor can be interpreted, if we 
consider the same integration as in Eq. ( |6?D in finite limits 0(Ao) = d\ 2 (dfi/d\ 2 )f 2 . 
The generated interaction, Eq. (|S6|) with 9(Ao)-factors, appears when high-energy modes 
are eliminated, since 9(Ao)-factor allows only momenta Q\ > Aq (and Q 2 > Aq). 
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The second term in Eq. (|65| ) is the instantaneous interaction, which comes from fixing 
of the light-front gauge || 

The sum of the dynamical and instantaneous terms is given by |5j 



6i , 6 2 \ SQ 2 Q l 



Q 2 Ql ^ l q+ 2 i g? Q 



2, 



where SQ 2 is given in Eq. (161]). In the light-front frame the energy transfers, Eq. 
read 



(x'k± — xk\ ) 2 + m 2 (x — x' y2 



Ql 



XX 1 

Ql = Ql\x-^x';{l-x)^(l-x') ? (70) 

that are always positive (> 0). 

We generalize the expression for an effective ee-interaction, Eq. (|57|) , to the case 
of QCD with a nonzero gluon mass. We simulate an effective interaction in QCD be- 
tween quark and antiquark as an exchange of non-perturbative gluon (gluon flux) with 
a nonzero effective mass (effective energy), which evolves with the cut-off A according to 
RG equation. The four momentum transfers read 



Q\{\) = Q? + /4p(A) 

+x';(l-a!)-»(l-a! / ) > 



Ql(X) = Q\{X)\ x ^ x '-ji-x)^(i-x>) , (71) 



where Q 2 are given in Eq. (|70|) ; and the nonperturbative effective mass, as obtained in 
the previous section Eq. Q5l|), is 

A 2 

VnpW = ~ <r(fi,v>) • (72) 

To reflect the phenomenological dependence of the effective gluon mass on the momentum 
we use the following parametrization 

Q 2 {\) = Q 2 -a\nE- 

Ql = Ql + ? 

a = cr(/2, u) , (73) 

which holds for Q 2 < jl 2 and Q 2 < A 2 . In fact it does not change the result for an effective 
gg-interaction what kind of parametrization to use,Eq. (|71f) with Eq. (f72|) or Eq. (|73|). 

In QCD we take into account the dependence of four-monetum transfers along quark 
and antiqurk lines on the cut-off. In full analogy with QED, the effective quark-antiquark 
interaction reads 

V q q = -Const{^Y)B^, (74) 
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where instead of coupling a em some constant term Const is introduced. At the end of 
calculations we fit Const and a to reproduce the correct coefficients for the short-range 
and long-range parts of gg-potential. Here B nv includes the dynamical and instantaneous 
gluon exchange. Eliminating by flow equations the quark-gluon coupling, where gluon has 
an effective cut-off dependent mass, one obtains the dynamical (generated) interaction 

B*? = D lv (q){Ii + h) , (75) 
where the integral term is given by 



df(Ql(\) /X 2 ' 
Ql(\) d\ 2 

with Q?(A) defined in Eq. (^J). 

In order to obtain the instantaneous interaction one should modify the QCD Hamil- 
tonian in the light-front gauge for the case of nonzero gluon mass. Instead, we use the 
same rules to calculate an effective gg-interaction as in the perturbative case for QED 
H . Applying flow equations to the light-front QED Hamiltonian, one can show, that the 
instantaneous propagator, fermion 7 + /2g + and gluon rfrf /q +2 , can be absorbed into 
the regular propagator, (tf+m) and D au (q) respectively, by replacing q, the momentum 
associated with the line, by 

? (i) = (V^ET- E Ti^ET-ET),^) , (77) 

\ in interm out in I 

in the numerator for those diagrams in which the fermion or gluon propagates are only 
over a single time interval. Here YLin iJ2out) denotes summation over all inital (final) par- 
ticles in the diagram, while Yl'mterm denotes summation over all particles in the intermedi- 
ate state other than the particle of interest. When energy is conserved, J2i n Q~ = J2 ou t Q~> 
we recover the rules of light-front perturbation theory as formulated by Brodsky and Lep- 
age ||. In Eq. (|77|) for index i — 1 the sign is plus, for i = 2 is minus. In order to absorb 
the instantaneous term Eq. (pSj ) and obtain an effective interaction Eq. (pD|), one should 
do the replacement in polarization sum of the dynamical term Eq. fl66|) 



E E^(?~ (i) )^> (78) 

i=l,2 i=l,2 



where are given in Eq. (|77|). 

Following this rule also for QCD, we obtain from the dynamical interaction Eq. ( |75f) 
an effective gg-interaction Eq. (|7l]), with 

SO 2 

£v = 9uu (h + h) + VnVv-^ (h - h) , (79) 
where the integral terms J, are defined in Eq. (|76|). 
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Similarity function is a function of the break Q 2 (X) / A 2 , where four momentum trans- 
fer Q 2 (X) introduces implicit dependence on the cut-off, Eq. ( |73f) . In this case the integral 
term, Eq. ([T6|), is reduced 



h 



d 



1 \ df( Zl ) 
X 2 ) dz x 



a 



(80) 



where z\ = Q\(\)/\ 2 and z 2 = Q\(\)/\ 2 . The first term, unit, describes perturbative 
one-gluon exchange (analog of photon exchange in QED); the second term arises from 
the dependence of effective gluon mass on the cut-off and defines the long-range part 
of the effective interaction, since it is more singular than the first term. Provided the 
properties for similarity function as in Eq. the integral factor Jj saturates by the 
values z\ ~ 1 and z 2 ~ 1, i.e. an effective range of integration is < A 2 < (A 2 , and A 2 ) 
where Aq ~ Q\ and A 2 ~ Q\. In this range of A the four-momentum transfers do not 
depend on the cut-off, i.e. from Eq. (|73|) Q 2 (A ) ~ Q\ and Q 2 (A ) ~ Q\. This enables 
to estimate the integral factor analytically. One can approximate the integral factor 
provided the following relation holds Q\ ~ Q 2 , which is considered below. 

We consider different similarity functions: exponential, Gaussian and sharp cut-offs. 
The corresponding integral factors read 



Qi + Ql \ + Ql 
Q\ + Qi\ Ql 




/ = exp(-Q 2 (A)/A 2 ), / 
/ = exp(-Q 4 (A)/A 4 ), /; 
f = 9(\ 2 -Q 2 (\)), 

where Q 2 are given in Eq. ([731). AVe define 

Q 2 = (Ql + Ql)/2 = Q 2 + v 2 
SQ 2 = (Q 2 1 -Ql)/2 = 5Q 2 , 

with Q 2 and 5Q 2 given by Eq. (|6l|) . The effective interaction between quark 
quark, Eq. ([79]), for the three choices of similarity function reads 



*1) 



(82) 
and anti- 



B 



fW 



B 



[IV 



(— + —\ i fa" Mt\ a _ 5 Q 4 _ 
9txv U 2 QV \Q 2 <?+ 2 J Q 2 Q 4 - SQ 4 

^U 2 QV U al Q 2 i +2 )Q 4 + 5& 



B 



fJIV 



9 " v W + Q\ 



Q 2 Q 2 Q 2 +5Q 



)) 



n +2 



a 



Q 2 





5Q 2 




Q 2 + 


5Q 2 



13) 
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where we defined Q 4 = (Q 2 ) 2 and 8Q 4 = (5Q 2 ) 2 . For 5Q 2 <C Q 2 the leading effective 
interaction in all three cases is the same and is given by the first term of Eq. (jS3|) 

V^-Cons t{ ^)(± + t) + (84) 

Gluon with an effective mass parameter jl mediates interaction between quark and anti- 
quark at the distances r ~ 1//2. 

We define the resulting gg-effective interaction in the limit of the gluon mass parame- 
ter fx — > 0. In this limit the average four-momentum transfer Q 2 , Eq. flS2|), and the string 
tension b = <j(fl,u), Eq. (|55"D, are given 

lim Q 2 = Q 2 

lim a = a . (85) 



From Eq. ( |8"4"D the leading behavior reads 



V q ,= -(ri»){C f a s ^+a^) , (86) 



which includes Coulomb and confining interactions. We show this explicitly below. Here 
Cf = T a T a = (N 2 — 1)/2N C . Here we restored the correct prefactors before the both 
terms, using the freedom to fit the overall constant, Const, and a term which is propor- 
tional to the scale u 2 , a ~ u 2 in Eq. (jpBj). 

We express the effective gg-interaction Eq. (|86|) in the instant frame. Instead of the 



light-front frame we use the instant parametrization p(p + ,k±) — > (p z ,k±), where the 
connection between the light-front x and z-component of momentum is given 



m 



In the instant frame the four momenta read 



Q = q -p zPz - 



2 




5Q 

where q = p — p' = (q z ,q±) is the three momentum transfer of the gluon. Here Mi and 
M.2 are the total energies of inital and final states, defined in Eq. ( |63D and the energy 
difference 5M 2 , defined in Eq. ([64]) , shows the 'off-shellness' of the process. In the instant 
frame one has 

M 2 = 4(p 2 + m 2 ) 

M 2 = A(p' 2 + m 2 ). (89) 



18 



that enters Eq. fl88|). As 5Q 2 — > the effective interaction has singularity at Q 2 — *> q 2 — > 0. 
The Fourier transform (with respect to q) reads 



which is the sum of the Coulomb and confining interactions. Though we were working 
in the light-front formalism the result for the leading effective interaction is rotational 
invariant. 



There are corrections O y^r) to the leading effective interaction, Eq. ( |S6D (or Eq. 

(|90|)). These corrections depend on the direction from which q approaches zero. For 
sufficiently smooth similarity functions f(z), as exponential and Gaussian cut-offs, the 
effective interaction does not contain collinear singularity (~ l/q + ). Thus the interaction 
becomes only singular if q approaches zero where it diverges like 1/q 2 (Coulomb singular- 
ity) and l/q A (confining singularity). However this is not true for the sharp cut-off, where 
liilv term diverges like l/q + - For a smooth cut-off the collinear singularity disappears 
(cancels completely) and only the rotational invariant part of the effective interaction 
survives in the limit q — > 0. 

5 Conclusions 

We suggested a possible scenery of confinement in the light-front QCD, basing on the 
method of flow equations. Flow equations operate in terms of 'physical' (dynamical) 
degrees of freedom, which are getting 'dressed' through the non-perturbative renormal- 
ization of the canonical QCD Hamiltonian. 

Integrating flow equation over the flow parameter in one gluon sector we obtained the 
gluon gap equation. It was solved, given an arbitrary mass parameter fi in the renormal- 
ization point, for the perturbative and non-perturbative gluon mass corrections. Perform- 
ing perturbative renormalization, the perturbative correction is absorbed by the second 
order mass counterterm. Severe collinear IR divergences, which arise in the gluon sector 
from the non-abelian gluon interactions, were regulated by introducing an additional 
cut-off u for the transverse momentum k± in IR region. This is equivalent to integration 
of gluon flow equation in the finite limits, from the bare cut-off A = Kjjv down to the 
hadron scale u, with u <C A uv . The result is the nonperturbative effective gluon mass. 

Eliminating the quark gluon coupling by flow equations, one obtains an effective 
interaction between quark and antiquark. In this approach the exchange with the dy- 
namical gluon mode between the probe quarks, where the effective gluon mass evolves 
with RG equations, gives rise to the effective gg-interaction which exibits Coulomb and 
confining singularities. The cut-off u, which regulates IR divergence, sets up a scale for 
the long-range part of interaction: it defines the string tension of confining interaction, 
a ~ u 2 . This suggests some relation between the zero modes of A + and confinement 
mechanism in the light-front formalism. 

For a smooth cut-off function the collinear singularity is canceled, and the leading 
effective interaction has rotationally invariant form. It is not true for the sharp cut-off. 




(90) 




19 



The ultimate goal of the study is to solve the chain of flow equations in different 
sectors selfconsistently. As was shown in this work, even an approximate solution of the 
gluon gap equation together with the flow equation for an effective interaction between 
probe quarks may provide an understanding of confinement. The next step is to include 
dynamical quark degrees of freedom, and to address in this formalism the problem of 
chiral symmetry breaking in QCD. 

This study shows, that in the light-front quantization it is possible to isolate the 
degrees of freedom that are responsible for the long-range properties of QCD, and ob- 
tain some insight into the non-perturbative QCD phenomena. Probably the light-front 
formalism is the most suitable frame to try solving selfconsistently the system of flow 
equations on computer. 
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A IR regularization via an effective gluon mass 

We consider correction to the effective gluon mass which arise from the non-abelian part. 
The IR singular behavior is regulated by the same mass parameter. From Eq. ( PB"D (the 
non-abelian part) one has 



A* a (A) 



/i 2 + 2g 2 C a I dx 



d 2 k ± 

JO 167T 3 

fi 2 (X)x(l-x) 



/(Qi(A)/A 2 



M 2 (A) 



k\ - x{l - x)n 2 {\) + (i 2 {\) k\ - x(l - x)n 2 (X) + /z 2 (A) 



1 + 



x- 



{i-xy 



(91) 



where similarity function, with momentum transfer Q\(X) given in Eq. (p4[) , regulates 
both UV divergences in transverse direction and IR divergences in longitudinal direction 



k 2 

Irani 



(X 2 + fi 2 (X))x(l-x)-fi 2 (X) 

v 2 W / _ , /Aa) 



A 2 + /i 2 (A) 



< x < 1 



A 2 + /i 2 (A) ' 



(92) 



Integration over k± gives 
/i 2 (A) 



9 2 C a 



X 2 x(l — x) 



., , dx U (A)x(l-x)(l + — )ln- 

7T 2 Jx ra ,; n V X 2 /i 2 (A)(l — X(l — X)) 



X 2 x(l — x) 



- M 2 (A)(l + ^)ln /x2(A)(1 _ x(1 _ x)) 
+ (l + ^)((A 2 + / , 2 (A))a;(l-a;)- yU 2 (A)) 



(93) 



where we have the symmetry with respect to interchange x and (1 — x). This may be 
simplified to 



/x 2 (A) 



fi 2 + 



9 2 Cg 

4tt 2 



/i 2 (A)ln^-(ln 



A 2 



Ai*(A) v /x 2 (A) 12' 



+ /AA)( 
1 



A 2 , A 2 
+ ln 



// 2 (A) ' "> 2 (A) 12 ; 1 12' 



dx(l + -=•)(! - x(l - x)) In 



x(l — x) 



1 — x(l — x) 



(94) 



We take into account the dependence of coupling constant on the cut-off, g(X), in Eq. 
) . Generally, the following terms contribute to the right hand side of Eq. (j94|) 



/ ^2 ^2 ^2 

H 2 {X) = fl 2 + g 2 (X) U 2 ( A) In— — (ci In ^^ + 02^^ + 03) 



V 2 (Xy /i 2 (A) > 2 (A) 



A 5 



V(A) 



(95) 
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where Cj and d { are some numerical constants. Following the same itterative procedure 
as outlined in the main text we substitute the leading order value for the effective mass 
into the right-hand side of Eq. fl95f), yU 2 (A) = p 2 . The next to leading order reads 

/ ^2 ^2 \2 ^2 \ 

/x 2 (A) =jl 2 + g 2 {\) fx 2 In — ( Cl In — + c 2 — + c 3 ) + fl 2 {c' 2 — + c' 3 ) . (96) 

V jJj j~-t j^t j. I J 

One may consider the following scenery, how the effective gluon mass fi appears in 
the theory. In our case the effective gluon mass plays the role of IR regulator. It is 
well known, that QCD (in chiral limit, m — > 0) initially does not have any scale, i.e. it 
is conformal invariant. When QCD is evolved from the bare UV cut-off to some lower 
energy scale, the canonical operators are changing with renormalization group equations 
in a way that the physical observables are expressed only through the renormalzation 
group invariant (cut-off independent) combinations. From Callan-Symanzik equation the 
RG invariant combination in the leading order of perturbation theory is 

A = Aexp bS))' (97 > 

so, dA/dX = 0. Here A is the running cut-off and A = Aqcd is the hadron scale, A A; 
b = fN c - |JV> is the one-loop Gell-Mann-Low coefficient in /^-function. The hadron 
scale, A, arises through the dimensional transmutation: one introduces the renormal- 
ization point - the parameter with the dimension of energy, A, in order to express the 
running coupling constant, g(X), through the coupling at renormalization point, g(A), 
which is dimensionless 0. 

We express all operators with dimension of energy through the hadron scale, since it 
is the only scale provided after renormalization. Therefore 

if ~ A 2 = AVp (-^) , (98) 

where in order to insure the boost invariance the cut-off is rescaled A — > X 2 /q + , and the 
same for A. The logarithmic term in Eq. (|96|) is given 

S 2 (A)lr4 = ^, (99) 
fl z b 



that reduces the Eq. (j96|) to 



A 



2 



At 2 (A) = ai jl 2 In — + a 2 jl 2 + a 3 X 2 + 0(g 2 (X)) , (100) 

where ai are numerical constants. To remove the leading cut-off dependence when A = 
A — » oo, we renormalize this equation by adding the mass counterterm 

m 2 CT = -a 3 A 2 . (101) 
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The nonperturbative effective mass is given 



VnpW — «iAi 2 ln — + a 2 /x 2 . (102) 

Indeed, this equation is of a nonperturbative kind, since it does not involve powers of the 
coupling constant. It was possible to get a nonperturbative result, because we regulated 
IR singularity by the hadron scale, which is given by a nonperturbative expression, Eq. 
fl97|): the argument of the exponent, 87r 2 /g 2 , can never be obtained in perturbation theory 
as an expansion in powers of coupling constant. 

In fact, we use a nonzero gluon mass jx as a parameter in our calculations, and at 
the end we take jl — > 0. However, from Eq. ( |102j ), when a mass parameter is removed, 
jl — > 0, an effective (nonperturbative) gluon mass vanishes (since jl is the only scale 
available through which li 2 np {\) is expressed). In the main text we introduce therefore 
an additional scale u to regulate IR divergences. 
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